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Flow rule, self-channelization and levees in unconfined granular flows. 
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Unconfined granular flows along an inclined plane are investigated experimentally. During a long 
transient, the flow gets confined by quasi-static banks but still spreads laterally towards a well 
defined asymptotic state following a non-trivial process. Far enough from the banks a scaling for 
the depth averaged velocity is obtained, which extends the one obtained for homogeneous steady 
flows. Close to jamming it exhibits a crossover towards a non local rheology. We show that the 
levees, commonly observed along the sides of the deposit upon interruption of the flow, disappear 
for long flow durations. We demonstrate that the morphology of the deposit builds up during the 
flow, in the form of an underlying static layer, which can be deduced from surface velocity profiles, 
by imposing the same flow rule everywhere in the flow. 
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Geophysical granular flows such as pyroclastic flows or 
debris avalanches self-channelize, forming a coulee sur- 
rounded by static banks, until they come to arrest and 
form a deposit [l|, Q with levee/ channel morphology. Re- 
cent laboratory experiments |3 involving short time in- 
terruption of a localized flow of dry granular material 
have reproduced such deposit morphology and under- 
lined the need for a deeper knowledge of the rheology. 
Flows down an incline provide a natural configuration for 
studying the rheometry of dense granular media. In the 
case of homogeneous flows, experiments |^ have shown 
that there is a minimum thickness /istop(^) below which 
no flow occurs and a maximum one /igtart {&) above which 
static layers spontaneously destabilise. Later it has been 
shown pl y, IM ly| that the depth averaged velocity uis 
related to the thickness h by the flow rule: 
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where g is the gravity, leaving unclear the dcpcndance 
of the non-dimensional parameters a and (3 on both the 
kind of grains and the covering of the incline. For a 
flow of glass beads down an incline covered with glued 
beads, a = and the flow rule is consistent with a local 
rheology |y, |3 ■ This rheology has been validated in var- 
ious flow configurations, including non-steady and non- 
uniform fiows as well as flows on erodible ground [a, |3 
and was recently extended to three dimensions |3| • How- 
ever, it does not verify u{hstop) = and hence does not 
describe the flow arrest, where non local effects are ex- 
pected to become significant [l^- Requiring this last 
condition imposes a = — /3, a situation actually reported 
in the case of unconfined flows |3- Hence, it is of pri- 
mary importance, both for practical and fundamental 
reasons to investigate further the situation where flow- 
ing and static regions coexist. 

In this letter, we study experimentally an unconfined 
flow down an incline. The flow indeed self-channelizes 




FIG. 1; (a) Sketch of the experimental setup, (b) Phase 
diagram in the plane {h, 9): the experimental data for /latop(^) 
(•) and hsta.rt{9) (o) with the best fits by eq. (solid line). 

within static banks which may evolve freely. Our main 
goals are to characterise the rheology and dynamics of the 
flow close to jamming, and its relation to the morphology 
of the deposit when the flow stops. 

Set-up - The experimental setup is sketched on 
Fig. ^a). A 60 X 300 cm^ plane, covered with sand pa- 
per of average roughness about 200 ^m, is inclined at an 
angle 6 ranging from 24° to 34°. The granular material 
consists of a slightly polydisperse mixture of spherical 
glass beads of diameter d = 350 ± 50 /im and density 
p = 2500 kg.m~ . Grains are released from a reser- 
voir located at the top of the plane by opening a gate of 
adjustable height and width. To produce homogeneous 
flows, the whole plane width is used, but for unconflncd 
flows the gate width is reduced to 5 cm and the width of 
the plane is always larger than that of the avalanche. For 
the present set up, /istart(^) and hstop{0) obey (Fig.[TlD): 
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with a = 3.2, 0-op - 22.5°, 0- rt = 23.0° 
and ^start = 35.6°. The mass flow rate Q is measured at 
the end of the plane, where the grains fall inside a reser- 
voir resting on a scale. During all experiments, some of 
which lasted up to 5 hours, Q fluctuated by less than 
2%. A 572 X 768 pixels camera positioned at the verti- 
cal of the plane is used to acquire images of the flow at 
25 Hz. The local flow thickness h is measured using the 
deviation of a laser sheet inclined at ^ 5° over the layer. 
The measured sensitivity (~ 60 ^m) is smaller than the 
grain size thanks to the averaging across the beam width 
(10mm) along x. About 10% of the grains arc dyed in 
black and the velocity of the surface grains, Us^ is mea- 
sured using a particle-imaging velocimetry algorithm. In 
all experiments the spanwise velocity component is found 
to be smaller than the resolution (2 mm/s). 

Transverse spreading & self- channelization- Upon re- 
lease of the grains, an avalanche front propagates down 
the plane at a constant velocity, leaving behind it a 
streamwise flow uniform in the x-direction (Fig. ^). 
Fig. [2Ia) displays thickness profiles at successive time 
steps, after the front has reached the downstream extrem- 
ity of the incline ~ typically after 100 s. It shows that 
the flow widens while the thickness proflles bend pro- 
gressively. The central flow thickness converges rapidly 
toward an asymptotic value Hoo{9, Q), whereas its width 
converges to an asymptotic value ^^00(6*, Q) with a much 
larger relaxation time (Fig.j^h)- When preparing a flow 
of width larger than Woo by increasing the flow rate, and 
returning to the initial flow rate, the width decreases back 
to Woo- This shows beyond any doubt the selection of 
an asymptotic steady state. Both Woo and Hoo are ob- 
served to increase with the flow rate Q. The aspect ratio 
-ffc5o/M^cx) turns out to be independent of Q and slightly 
decreases with 9 (Fig.Et). 

Fig- 13d) shows the transverse structure of the flow: an 
inner flowing region (u ^ 0) flanked by two static banks 
{u = 0) of width approximately 5 mm. These banks 
forms due to a strong increase of friction on the side of 
the layer, as the thickness h vanishes and thus /i tends 
to Ho = tan 00 (eq. @). Once the banks are formed, 
their external sides arc so steep that the free surface an- 
gle tan$ = (tan^(0) -I- {dyh)'^^/'^ becomes significantly 
larger than the plane angle (Fig. [2t) allowing trans- 
verse displacements. Indeed, the banks are outside the 
metastable band [/istop, ^start] during the spreading phase 
and converge toward it in the asymptotic state. The flow 
is hence divided into a central flow region uniform in the 
y-direction - the ratio of the transverse to the vertical 
diffusion of momentum h? /W^ « 10~^ - and the banks 
dominated by three-dimensional effects. In both regions, 
incrtial terms Uydux/dy w hul/W"^ « 10~'^5 are small. 

Flow rule - The flow rule is flrst determined using ho- 
mogeneous steady flows covering the total width of the 
plane (Fig.|2t, A). For h/hstop > 2 it is well described by 
eq. (Q) with a — and /3 = 0.134 assuming u/ug = 3/5 




3- 







*y(«n>^'' 



FIG. 2: (a) Time evolution of the height profiles h{y, t) and 
(b) evolution of the central height H and width W of the 
flow for Q — 2h g.s^^, 6 — 25°. (c) Asymptotic aspect ratio 
of the flow Hoo /Woo as a function of the slope angle, for flow 
rates Q ranging from 3 g/s (A) to 30 g/s (a). (d) Velocity 
Us{y) (solid line) and thickness h{y) (dashed line) proflles for 
e = 25°, Q = 25 g.s"^ and t = 150 s. (e) Thickness vs. local 
slope $ inside the quasi-static banks (o) : f < 2000 s, (•): 
asymptotic state. The solid lines show ftstop and ftstart. 

(see below). For h/hstop < 2, there is a systematic devi- 
ation towards lower Us/{gh)^''^ values. However smaller 
thicknesses can not be further investigated in homoge- 
neous flows. On the contrary self-channelized flows are 
naturally very thin and give access precisely to that range 
of thickness. Fig. |3t (•) shows data collected near the 
centerline of the flow, where the flow has been shown 
above to be uniform in the transverse direction. Plot- 




FIG. 3: (a) Flow rule Us/ighf^ vs h/h^top- (A) data for 
steady homogeneous flows, (•) data obtained in the central - 
uniform - part of unconfined flows. Each point is an average 
over different realizations and the corresponding statistical 
error bars are reported. The solid (resp. dashed) line is the 
best fit by eq. Q (resp. eq. Q with a = 0). (b) Average 
to surface velocity ratio u/ua vs. h/hstop for steady homo- 
geneous flows, (c) Scaling of the asymptotic height with the 
flow rate: Qh = pg^ HJ, {Hac/hstop — 1) as a function of 
Q, for e ranging from 25° (A) to 28° (A). 

ting Ual \fgh as a function of h/hstop, one observes again 
a collapse of data (for different 9, Q and t) on a single 
curve which satisfies u{hstop) — and can be described 
at the first order in h/hgtop by: 
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with j3 ~ 0.219 assuming now u/us = 1/2 (sec below). 
Note that both in the confined and unconfined cases, we 
have checked using dyed grains that the fiow involves 
the whole layer thickness. Most importantly, the data 
obtained in both cases coincide on their common range 
of thickness and provide - to the best of our knowledge - 
the first experimental determination of the flow rule over 
the whole range of thicknesses, independently of the flow 
configuration. 

We now briefly discuss its relation to the locality of 
the rheology. For a large enough thickness, the flow rule 
obeys eq. (^ with a = and is consistent with a lo- 
cal rheology p : when the shear stress r depends on 
the shear rate 7 at the considered point only, it follows 
from dimensional analysis - assuming a small transverse 
shear stress as suggested by the flatness of the trans- 
verse velocity profiles - that the rheology can be ex- 
pressed under the non-dimensional form t/P = fii{I) 
with / = 'yd/{P/py^^, where d is the grain diameter, 
P the pressure and p the density. Integrating this last 
relation leads to a scaling law for the depth averaged 



velocity u(h,8) — A{9)h^/'^, which is indeed satisfied 
by eq 1^, when a = 0. For a small thickness the 
flow rule eq. ((SJ clearly displays a violation of the lo- 
cal rheology close to jamming, which is confirmed by the 
measurement of the ratio of the depth averaged veloc- 
ity u - measured with the avalanche front velocity in 
the homogeneous flow case - to the surface velocity Us 
(Fig. Ob)- For large h/hstop values, u/us = 3/5 which 
is consistent with Bagnold-like profile deriving from a 
local rheology [a|. Close to jamming transition, u/ug de- 
creases to 1/2 in agreement with numerical findings [5J 
of a transition toward a linear velocity profile. Note that 
the above analysis suggests that P/hstop{0) shall not de- 
pend on the covering of the incline whereas /3/ft.stop(^) 
and of course /istop could. Finally from conservation 
of mass, one expects a linear relationship between Q 

and HooWoou{Hao) oc pg'^/'^HoL (^oo/^stop - 1), using 
eq. (PJ. This relation is closely verified experimentally 
(Fig. Ot) and suggests to describe the flow in a shallow 
water approximation. However, let us mention that nei- 
ther the asymptotic state (its bended shape but also the 
selection of the aspect ratio H^/Woo) nor the temporal 
scalings of W{t) and H{t) can be simply predicted fol- 
lowing such an approach. Hence we shall postpone this 
analysis as well as a full three-dimensional analysis to a 
forthcoming paper. Note that the complete set of shallow 
water equations is solved in |ll| . 

Banks - In the previous section, we investigated the 
flow rule near the centerline of the flow. No such scal- 
ing is a priori expected close to the static banks. Indeed 
as seen on Fig. 0Jb), when acquiring data moving away 
from the centerline, no scaling is observed and Us/{gh)^''^ 
is systematically underestimated. A possible explanation 
of such a deviation from the otherwise scaling is the ex- 
tension of the static zone below the flowing grains. In 
order to verify such an hypothesis, we have performed 
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FIG. 4: (a) Static layer thickness Z reconstructed from sur- 
face velocity measurements (solid line) and directly measured 
with a brush (o) - dashed line : total height, (b) Ua/{gh)^'^ 
vs h/hstop, data obtained across the whole flow width: there 
is a progressive shift from the fiow rule when moving away 
from the centerline (arrow). 
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FIG. 5: (a) Thickness profiles of the deposit /idcp(y) obtained 
after different flow durations for Q = 25 g.s^^,9 = 26°. 
(b) Comparison between measured (solid line) and predicted 
(dotted line) deposit profiles, (c) Picture of levees observed 
at the borders of a deposit. 

specific experiments. 

One first conducts an experimental run using white 
grains. After a few minutes, the flow is stopped and a 
deposit forms. The white grains of a small slice across 
the deposit are removed and replaced by black grains 
of the same material. The flow is then started again 
at the same flow rate. At the surface, the black grains 
are washed out by the incoming white grains except in 
the static borders. After a few minutes, the flow is in- 
terrupted again. Using a brush, the white grains are 
removed very cautiously layer after layer in the region 
of the slice. Close to the centre, the deposit is exclu- 
sively composed of the new white grains, proving that 
the flow involves the whole thickness. However, close to 
the banks, a layer of black grains remained trapped, in- 
dicating the presence of a static [l^l layer of thickness 
Z{y,t) (denoted by o in Fig.^) below the flowing one. 
The interface between black and white grains turns out 
to be very sharp (of the order of a single grain diameter) . 
Using Us, one can construct the thickness R < h which 
would flow if the flow rule obtained above were applicable 
inside the flowing layer. The reconstructed static/flowing 
interface Z = h — R is shown on Fig. EJa) together with 
the measured one. The fairly good agreement suggests 
that, to the flrst order, mobile grains flow above static 
grains just like they flow above the rough plane. At the 
same order, it tells us that the flow rule is actually ob- 
served everywhere in the flow when applied to the flowing 
thickness R instead of the total one h. 

Levees - Turn now to the morphology of the deposit 
which forms when the flow stops. It has been reported J3 
that, under certain conditions of inclination and flow 
rate, the deposit formed upon interruption of the flow 
exhibits a levee/channel morphology similar to those ob- 
served on pyroclastic flow deposits. In the present study. 



we have observed that for larger times, the flow actually 
keeps on widening and converges only very slowly to- 
ward its asymptotic state. Accordingly the shape of the 
deposit strongly depends on the flow duration t (Fig.|SJi). 
For small t, deposits are composed of a central flat zone of 
thickness /istop bordered by two levees of thickness larger 
than /istop (Fig. 13;) as previously reported j^]. When t 
increases, the levee thickness decreases until it vanishes 
at very long time, so that the deposit corresponding to 
the asymptotic state is indeed flat. Levees result from the 
combination between lateral static zones on each border 
of the flow and the drainage of the central part of the 
flow after the supply stops [3j. However, a clear picture 
is still lacking concerning the junction between a central 
flat zone of thickness /igtop and a levee of thickness larger 
than ft-stop- The flow rule that was obtained here provides 
a very simple scenario. The flow stops when R = hgtop, 
and one expects the deposit to form by superimposing a 
layer of thickness ft-stop to the static layer. Accordingly 
the deposit thickness should be /idep = Z + /igtop wher- 
ever Us ^ and /idcp = Z in the quasi-static banks, where 
Us = 0. Fig|S{b) provides the experimental evidence that 
such a simple scenario indeed holds: the predicted and 
the measured deposit profiles qualitatively match. 

To conclude, by investigating unconfined granular 
flows down an incline, we have shown that they obey 
a flow rule which reveals a crossover towards a non lo- 
cal rheology close to the jamming transition. This flow 
rule accounts for the morphology of the deposit, which 
actually builds up during the flow, in the form of an un- 
derlying static layer. This underlines the importance of 
addressing erosion-deposition mechanisms issues to com- 
plete a full description of geophysical flows. 
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